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ABS'IRACT 


The water intake devices of Kanpur city’s water works 
located at Bhairoghat on the right bank of River Ganga had 
adequate water supply till the mi d-twent i e th century. However, 

gradual migration of the river towards the left bank since the 
1950*s produced acute water crisis at water intake structures 
during the lean period of flow. This necessitated digging and 
maintaining a 6 Km long uniined canal to bring water from the left 
bank to Bhairoghat- Moreover, shifting of the water course 
towards the left bank caused impingement of flow and scouring at 
the earthen embankment of Kanpur -Lucknow road and railway bridge 
near Champapur. Both these problems may be rectified if the river 
course near Kanpur shifts towards the right bank. 

A hydraulic model study was attempted here to train the 
river to flow near the right bank using permeable spurs for river 
training. A distorted scale model -- horizontal 1 in 9090 and 
vertical 1 in 130 -- with mobile alluvial bed was designed. Ganga 
river’s flow data of 28 years, bed material data collected on 
site, and river planform and bed level data were used to simulate 
the river dynamics. Copper rods Anvti diameter, 4 to 5 diameters 
apart, were spaced in sequence to simulate proposed spur arrays in 
the river. Average and maximum annual flow hydrographs of 28 
years were used to simulate the flow conditions in the river, 
using various types of spur array patterns. It was observed that 
a double array of permeable spurs near the left bank, and at acute 
angle to the bank upstream of Bhairoghat, was extremely effective 
in altering the river thalweg profile, shifting the river course 
towards the right bank and reducing scour at embankment. The 
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study indicates very positive results using permeable spurs for 
river training, but suggests the need for a more elaborate and 
detailed model study. 
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CHAPTER I 


INIRODUCnON 


1.1 

GANGA 

River 

RIVER PROBLEMS 

Ganga meanders 

NEAR KANPUR 

or braids in 

most places 

i n 

the 

plains 

of Uttar 

Pradesh. Near 

Kanpur city. 

the river has 

a 

wi de 


section but the thalweg profile, marking the water course during 
the lean summer months, meanders far away from the city towards 
Unnao. This was not always so; till 1956, the flow path was 
Kanpur side, and the water intake structures near Bhairoghat could 
draw enough water for the city’s major domestic and industrial 
water needs. However, the sun'mer water course having shifted 
towards Unnao, water is now brought to the intake structures by a 
6 Km long cl»anrH'‘l. The regular maintenance and deal 1 ting 
operations of this channel being costly, the city’s engineers and 
planners have long sought a viable alternative to this 
arrangement. Figure 1,1 gives the location of River Ganga near 
Kanpur. 

A second related problem is the severe scour at 
embankment near Champapur built for the road and railway bridges 
connecting Lucknow and Kanpur. The bridges, located about 5 Km 
downstream of the city’s water works are on the extreme right of 
the flood plain, whereas the thalweg profile i rrwed i at e 1 y upstream 
of the bridges is on extreme left as can be seen in Figure 1.2 
which gives the pianform details of River Ganga near Kanpur. 
Thus the flow turns sharply skirting the guide banks* before 
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turning left again, almost at right angles, to flow under the 
bridges. These sharp changes in flow direction cause excessive 
scour near the nose of guide bank. During high flood the severe 
scour causes danger to the embankment. These problems prompted 
the model study of River Ganga near Kanpur. 

l.S OBJECT OF STUDY 

The aim of this work is two-fold. First, is to make 
available sufficient flow during the lean period at Kanpur city’s 
water intake structures; second, to protect the guide bank of 
bridge and embankment from the severe scour. Both these objects 
were hoped to be achieved by a simple and relatively inexpensive 
river training method. 

A laboratory study was pl«aru>ed t.o actueve the above 
Dbjects. A distorted mobile bed model from upstream of Bithur to 
the road and railway bridges was planned. The use of permeable 
spurs was contemplated to divert the lean flow towards Kanpur 
3ank • This was expected to ensure lean period flow at or near the 
aater intake structures and also to avoid the severity of scour 
tear railway bridge guide bank and embankment both during lean and 
'lood flow periods. 

.3 PRESENTATION OF THE REPORT 

Ihe thesis is presented in six chapters. The first 
hapter introduces the problem and the thrust of the present work, 
he second chapter briefly describes the relevant concepts and 
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knowledge of river modelling based on literature. In Chapter 5 a 
comprehensive description of the prototype problem, historical 
information, morphological and hydrological features are 
presented. The fourth chapter documents the prototype data used 
and their sources, and the design of experiments. In Chapter 5, 
the experiments have been described in detail, and the results 
analyzed. The important conclusions and limitations of this work, 
and suggestions for further work are listed in the concluding 


chapter (Chapter 6). 



CHAPTER II 


LITERATURE SURVEY 


2*1 INTRODUCTION 

A river is a natural water course carved on land over a 
long period of time* and conveying water from upstream sources 
(reservoirs, icemelt, etc.) or from precipitation in a catchment 
(surface run-off or subsurface flow) to a downstream surface water 
body. The diversity of factors and the complexity of their 
interactions that govern a river’s contemporary status -- the 
morphology of the land, the extent and distribution of its water 
sources, the history of its development, the effect of manmade 
developments, the mechanics of water movement, the sediment 
carried into it by overland flow and strong winds, ambient 
conditions such as temperature, and even rotation of the earth on 
its axis makes it a most complex process for scientific study. 
While a river may be thousands of years old and a thousand 
kilometers long, the hydraulic engineer’s interests and knowledge 
are usually limited to much shorter life-spans and lengths of 
rivers. Thus the hydraulic river researcher can delineate several 
important variables that characterize its regime, such as water 
and sediment discharges, valley slope, roughness, its 
three-d i mens i ona I geometry (namely planform, and cross sectional 
dimensions). While conceptually the flow processes in a river may 
be well understood, the actual quantitative relations between 
different variables are difficult to determine. In practice, a 
large number of empirical or semi emp i r i ca 1 relations (such as 
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Manning’s resistance equation, sediment transport relations, river 
regime formulae) are used to define the river regime and to 
predict changes in river behaviour [Leopold and Wolman. 1964. 
Jansen, et al . 1985; Knighton, 1984], Consequently, mathematical 
formulations are not often adequate, and physical (scale) models 
may be used to predict changes in river behaviour, particularly 
before attempting to modulate or control actual river behaviour. 

2*2 RIVER M0RPI10LCX9Y 

Rivers are characterized by both time varying and 
spatially varying quantities. The water and sediment discharges 
vary seasonally ar'id annually. The river course, on the other 
hand, is usually non-un i f orm. Planform studies on rivers have 
delineated three major patterns, --- straight, meandering and 
braided, as illustrated in Figure 2.1. Straight rivers are found 
mostly on steep slope, but hydraulically such rivers may not be 
uniform due to the formation of step-pool systems or other 
non-uniformities in cross section. 

Meandering rivers are very common on alluvial plains. 
When once a river deviates from its axial path, a curvature is 
developed (either due to its own characteristics or due to the 
impressed external forces) and the process moves down-stream by 
building up shoals on the convex side by means of secondary 
currents. The formation of shoals on the convex side results in 
further shifting of the outer bank by eroding on the concave side. 
Formation of successive bends of reverse order may lead to the 
formation of a complete S-curve called meander. When consecutive 
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curves of reverse order connected with short straight reaches, 
called crossin gs» ar^e developed in a river reach then the river is 
stated to be a meandering river. Some useful ways of 

characterizing a meandering river are given below. 

Thalweg t It is the line of maximum depth of stream along 

which the flow moves in lean period also, 

Meandor longih C ) i It is the axial length of one meander, 

i.e., the tangential distance between the corresponding points of 
a meander . 

Meander belt or meander width t It is the distance between 

the outer edges of clockwise and anticlockwise loops of a meander. 
Meander Ratio ^ ^ * It is the ratio of the meander width to 
the me a n d e r I e n g t h . 

Sinuosity : It is defined as the ratio of thalweg length to the 

valley 1 ength . 

Tortuosity t li-i M x 100 

va 1 1 ey 1 ength 

Braided rivers or streams may also be observed on 
alluvial plains, and they are corrmon near tidal outlets of rivers. 

A braided river can be defined as one which flows in a number of 
channels around alluvial islands. Figure 2.1 shows a typical 
braided reach of an alluvial stream. 

Laboratory studies on river patterns have been conducted 
extensively in this century, particularly since Friedkin's study 
on river meandering C1945]. Field observations also have led to 
considerable insight and knowledge of river morphology in relation 






SKETCH FOR MEAND E RINO AND BRAIDED 


Figure 2.1 


Definition 

STREAMS 
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to some dominant variables (Leopold et al., 1964, Knighton, 1984]. 
However, neither is the mechanism of morphological adjustments of 
a river fully understood, nor is there a reliable method to 
predict morphological patterns and changes accurately, 

2.3 RIVER MODELLING 

2# 3# 1 Scale models and mathematical models t 

There are two main types of mobile bed models 
mathematical model and scale model (or physical model). Analog 
models are of little use since they can not mimic the complex 
mechanism of river adjustments. 

(i> Mathematical model : Certain problems in river engineering 
can be conveniently investigated by means of mathematical 
models. It is a model using mathematical relationships to 
represent the prototype. A fundamental pre-requisite of 
mathematical model is a satisfactory quantitative 
mathematical description of the physical processes involved. 
This description very often consists of a system of ordinary 
partial differential equations together with suitable 

boundary conditions and other data. 

<ii> Physical model : A physical model is a precision device used 
in order to predict the behaviour of a physical phenomenon. 
Unlike analog model, it uses the physical properties and 
behaviour of modelling materials to represent the 

prototype [ ASCE Task Corwii 1 1 ee , 1982 ]. The principal 

difference between mathematical modelling and physical 

mode 1 1 i ng cons i s t B in the fact that a mathematical model 
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requires the formulation of equations which describe the 
phenomenon or flow field whereas it is sufficient for the 
hydraulic model to identify the acting forces and from these 
to formulate similarity parameters. 

The relative merits of the above two types of models may 
be described as below. 

(i) From the numerical point of view three dimensional time 
dependent computation will be out of the question for some 
time, at least for practical purposes. But solution of two 
dimensional time dependent problems are within the capability 
of existing computers. Physical models can be used to solve, 
both two and three dimensional problems, but their 

I 

applicability is selective depending on the hydraulic 
cond i t i ons . 

(ii) From physical point of view the process of turbulence is not 
yet fully understood. This means that the formulation of 
momentum transfer or mass transfer by turbulence contains 
empirical elements which limit range of application. 

If one is faced with the decision to solve a problem 
either by means of a hydraulic model or . by means of numerical 
model, one has to consider a variety of aspects as criteria in the 
decision process. The consideration of principal limitations may 
exclude a priori the one or the other type of model for certain 
problems, and which degree of accuracy or resolution is required 
from the model. The most important role in the decision making is 
played undoubtedly by the limiting factors of either type of 
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mode i . 


Table Z.l gives a comparative chart of the factor 
limiting the usage of hydraulic (scale) and numerical 
(mathematical) models [ Kobus . 1980 ]. 

Table 2*1 Limiting factors for live-bod models 


Hydrau 1 i c mode 1 

Numerical model 

Pr i nc i pa 1 

L i mi tat i ons 

Model size (laboratory) 

Discharge (pumping capacity) 

Energy head (pumping capacity) 

Model laws 

Storage Capacity 

Computational speed 

Incomplete set of equations 

Turbulence hypothesis 

Prac t i ca 1 1 

Limitations 

Minimum model scale 

j 

In simplified set of equations : 

-accuracy of assumed relationships 

-availability of coefficients 

Model size (upper limitations) 

-space and time resolution 

(lower limitation) 

Measuring method and data 

coll ec t i on 

Numerical stability and 

convergence of the solution 

scheme 

Availability of boundary and 

initial cond i t i ons 

Availability of boundary and 

initial cond i t i ons 


2# 3. 2 Physical models 

A model in its widest sense is a simplified 
representat i on of a subject, state or event. Models may be 
similar or d i s-s imi lar . In similar models, all model parameters 
exhibit a certain relationship to the corresponding parameters in 
nature or prototype, which is determined by one or several model 
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scales. In d i s-s imi lar models the above requirements is not or 
only partially satisfied. 

Hydraulic models are small scale reproductions of nature 
in the laboratory. In hydraulic model tes t i ng , usual ly model fluid 
used is water due to its being easily available, cheap and simply 
replaceable. However, the similarity requirements do not prevent 
the use of other fluids and may be advantageous in some cases. 
For example, similarity of sand transport in waterhas been 
achieved in small scale models by means of the transport of coal 
dust in gl ycer in. 

2# 3* 3 Similarity mechanics 

The model and prototype should have similar properties. 
And this similarity between model and its prototype is called 
similitude. For absolute similitude between a model and the 
prototype, the following types of similarities should exist. 

i> Geometric similarity : This is the similarity of form which 

exists between model and prototype, if the ratio of 

corresponding linear (cartesian) dimensions of the model and 

prototype is constant. If the geometrical length is in 

nature and corresponding length in model is then length 

scale number I = I /I . 

r n m 

ii) Kinematic similarity : This is the similarity of motion 
which exists when the ratio of corresponding kinematic 
quantities at corresponding points of the model and prototype 
are the same i.e., ve 1 oc i ty , acce 1 erat i on , etc. 
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Velocity scale number V = V /V 

r n m 

Acceleration scale r'iumber a = a /a 

r n m 

iii> Dynamic similarity : This is the similarity of masses and 
forces which exists between model and prototype when 

(a) The ratios of masses of . correspond i ng fluid particles in 
motion are the same. 

<h) Thn ration of forrnn of^ ro r r n npond i ng fluid particln© 
are the same. Dynamic similarity includes geometric and 
k i n ema tic similarities. 

If Fn is the force acting in nature and Fm is the 
corresponding force in model, then force scale number F^ is 

F = F /F 
r n m 

Dimensional analysis allows the formulation of criterion 
for dynamic similarity, since all relationships derived by 
dimensional analysis are independent of the absolute scale, they 
must be applicable for both small scale model and prototype 
dimensions. 

A dimensional analysis leads to the conventional and 
well-known fluid mechanics characteristic numbers and model laws. 
By considering a fluid element of density p with a referehce 
length I and reference velocity v under the influence of viscous 
forces, one obtains by dimensional analysis a fluid mechanics 
characteristics number as a ratio of inertial and viscous forces 


Re 


Inertial reaction 
viscuous force 


pvl 


< 2.1 > 
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Re is known as Reynolds number. Very small Re characterize by 
definition flows in which viscous forces dominate and inertial 
reactions are negligible, as for instance in Darcy flow; flows in 
which the viscous forces are negligibly small in comparison to the 
inertial reaction, as for instance in fully turbulent pipe or 
channel flows, have high Re. 

Cor respond i ng dimensional considerations for a fluid 
element with inertial reaction under the influence of gravity 
forces yields another characteristic number, Fr. 


1 h e f o r' c e r a t i o f r 


i no r t i a I reaction 
gfavTly f orFci" 


V 

■/gh 


<2 . Z) 


Fr is called Froude number, which plays a dominant role in 
hydraulic modelling, as for instance in flows with a free surface. 
Negligible influence of gravity in comparison to inertial reaction 
corresponds to very large Froude ■ numbers, whereas very small 
Froude numbers correspond to an overwhelming influence of gravity 
forces . 


2.3.4 Model laws 

Reynolds model law : In flows with significant viscous effects. Re 
has to be kept the same in model and prototype. Apart from 
geometric similarity, this requires in addition satisfaction of 
the Reno.lds model law 


Re 

n 

Ra“ 


m 




1 
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This law is satisfied if the velocity in the model is 
chosen such that their results 


V 

n 

m 


I 

m 


7“ 

n 




or V 

r 



r 


<2. 5) 


If the model and prototype fluid are same, as in many laboratory 
experiments, the velocity number becomes 

This implies that in a small scale model the resulting 
velocities must be larger than in the prototype. 

Froude model law i For grav i ty -dr i ven flows geometric similarity 
and Froude number must be preserved in model and prototype. Hence 
Fr = V /Vq h =1. Sinceg =1, 

1"'* I' I' X*' 


V 

r 



(2.4) 


2* 3# S Fixed bed models 

''Fixed bed’ refers to a bed whose roughness 
characteristics do not change due to the flow. In general models 
whose beds and banks can not be eroded are called fixed bed 
models. Such models are useful in the investigation of problems 
covering long regions of a river in which variations in bed 
configuration are not considered. Such models are used in the 
following studies : 

i> To study changes caused by placing obstacles to flow like 
dams, piers etc. 

ii> Effect on navigation conditions and backwater profiles during 


f 1 oods . 
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Pi started and undi storie d mode 1 s : To model water movement, it is 

necessary to ensure Froude number similarity. In an undistorted 

model the vertical and horizontal scales are identical, i.e. I = 

r 

h , and according to Froude relation 
r 


V - h or I and Q 

r r r r 


= I 


5 /: 


The model is perfectly defined by the only choice of the 
geometric scale number Then it is necessary to verify that 

frictional forces are scaled in the same manner as the inertial 
reactions. It means three conditions have to be fulfilled. 

i) The Froude number must be the same in model and prototype. 

ii) The roughness of the model must be correct. 
iii> The flow in the model must be turbulent. 

The difficulties conrmonly observed are : 

i) Problem of measurements because depths and velocity will be 
sma 1 1 . 

ii) The Reno Ids number in the model is not high enough to 
guarantee a turbulent flow, unless model fluid is different. 

iii) The roughness of the model may not conform to prototype 
roughness . 

To overcome the above difficulties in a scaled down 
model a solution consists in distorting the model by exaggerating 
the model depth. In these distorted models it is necessary to 
distort the slope in order to overcome the high resistance and 
also to provide a suf f i c i en t I y high value of the Reynolds number 
to ensure turbulent flow conditions. 
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as 


Distortion of slope can be estimated from Manning’s law 


and 






< 2 . 5 > 


( 2 . 6 ) 


2*3.6 Mobile bed models 

These are models of channels and streams, meant for 
investigating problems involving erosion, transportation and 
deposition of material of the channel bed. These models are very 
difficult to design. It is known that bed load is moved due to 
the tractive force exerted by the stream. When geometrically 
similar models are made, the tractive force available is usually 
insufficient to produce corr espond i ng bed material movement when 
sand is adopted as the bed material. Hence in order to provide 
adequate tractive force mobile bed models are generally made 
distorted so as to provide steeper slopes. It is possible to bring 
down the value of the tractive force and, hence the extent of 
model distortion, if lighter materials of specific gravity greater 
than unity are adopted, eg. bakelite powder, pumice, coal dust, 
stone dust, wax balls, glass beads, plastic grains, etc. 

It is also important that the scale selected is such 
that it provides a sufficiently high value of Renolds number to 
ensure turbulent flow, even at minimum discharge condition. 
Distorted mobile bed models : Models of rivers and harbours are 


usually constructed by different scales for horizontal 


and 
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vertical dimensions to remove possible surface tensile effects and 
to maintain turbulent flow condition as exists in the prototype. 
Hence vertical scale is exaggerated resulting in a distorted 
model. In principle it is hydraulic similitude and not geometric 
similitude, which is the main governing factor in model design. 
Hence it may be necessary to adopt distorted models. The 
following types of distortion may be adopted. 

i> Geometrical distortion : This is a distortion introduced by 
adopting different scales for horizontal and vertical 
dimensions. 

ii) Configuration distortion : In this case the bed slope of the 

model is increased; otherwise the model is geometrically 
s imi lar . 

iii) Hydraulic distortion : In this type of distortion some 

hydraulic quantity, say velocity or discharge, may be 
changed . 

iv) Material distortion : This type of distortion involves the 
use of material different from those in the prototype. 
Surface roughness or the medium in which the model works may 
be changed. 

Mer i t s of d i s tort ed mode I s : 

i> Higher Re in the model (especially to avoid laminar flow in 
f lood plains ) . 

ii) Improvement of relative accuracy of measurements. 

iii) Shortening of model test durations, 

iv) Decreasing of the importance of water losses in the model. 
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Demer i t s o f distorted f no dels : 

i) Due to unequal horizontal and vertical scales the pressure 
and velocity distributions are not truly reproduced in the 
mode 1 . 

ii> The wave pattern in the model will be different from that in 
the prototype due to depth distortion, 
iii) Slopes, bends and earth cuts are not truly reproduced. 

2# 3* 7 Ri ver model s 

For a river model with mobile bed the similarity 
requirements for .sediment transport must be satisfied in addition, 
and the simulation of transport process requires furthermore the 
.proper reproduction of density currents, turbulent dispersion 
process and reaction and exchange processes. 

Water flowing in rivers transports solid material of 
various kinds and rates. Whereas the upper river reaches contain 
mostly coarse material, which at the larger slope is transported 
as bed load, the lower reaches contain usually finer materials, 
the motion of which is mostly in suspension. The sediment 
transport, which varies with discharge conditions, causes changes 
in bedf orms , which can affect changes of the entire river course 
in the untrained river. 

A great number of different i nves t i gat i ons have been 
devoted to the solution of sediment transport problems. Whereas 
in the former times mainly empirical approaches were attempted, 
theoretical i nves t i gat i ons have led in recent years to an. 
essential expansion of our knowledge of sediment motion. Although 



the sediment transport equation allows a rough estimate, it is 
usually not possible to calculate accurately the mutual 
interaction between river bedforms and discharge distribution. 
Therefore, the solution of difficult problems of river hydraulics 
requires often the use of models with mobile bed. Whereas in 
model tests with fixed bed the Froude model law is sufficient, 
model investigations with mobile bed must satisfy additional 
conditions with regard to the horizontal and vertical length 
scales and with regard to the bed materials. Procedure for 
obtaining dyruimic (mobile bed) river model scales using 
dimensional analysis have been outlined and further discussed in 
the appendix. 

2.4 TRAINING WORKS AS RIVER AID 

River engineering works may vary greatly in size and in 
their effect on a river’s behaviour. They may provide a strictly 
local improvement which hardly changes in upstream or downstream 
river reaches; or they may be designed to completely alter a river 
regime, thus affecting the river over almost its entire length. 
Whatever their scale may be, however, river improvement schemes 
always need careful consideration and a scientific as well as 
practical approach. 

River training, in its wider aspects, covers all those 
engineering works which are constructed on a river, so as to guide 
and confine the flow to the river channel, and to control and 
regulate the river bed configuration, thus ensuring safe and 
effective disposal of floods and sediment loads. Stabilizing and 
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raining the river along a certain alignment with a suitable water 
ay in tiie first arid foremost aim of river training. The ultimate 
im of river training is to achieve permanent stability of the 
iver. Here, stability of a river rrreans that the river attains an 
quilibrium state, and no significant change occurs in its 
lignment, slope, regime etc. 

River training works commonly involve groynes and spurs, 
here are two types of spurs used in the field, namely Ca> 
mpermeabie spurs, and (b> permeable spurs. 

a) Impermeable spurs : are also called solid spurs or 

embankment spurs. These spurs may be rock fill embankments, 
armoured with stone pitching, concrete blocks etc. These 

spurs are called ’M mperemeabl e” because they do not allow any 
significant flow through them. 

[b> Permeable spurs : Permeable spurs do permit restricted flow 
through them. These spurs are more or less temporary 

structures and are susceptible to damage by floating debris 
etc. Permeable spurs simply obstruct the flow and reduce its 
velocity, causing silt deposition in the vicinity. They are, 
therefore, best suited for rivers carrying huge sediment 
load. Permeable spurs do not change the flow abruptly as is 
done by impermeable spurs, and hence, intense and serious 
eddies and scour holes are not developed. 

In this study permeable spurs are used as river training 
aid to control alignment and bank scouring. 



CHAPTER 3 


FEATURES OF THE PROTOTYPE PROBLEM 

3-1 HISTORY OF COURSE OF GANGA RIVER SINCE 1910 NEAR KANPUR 

The meander pattern of Ganga river in the twentieth 
century has changed considerably as shown by records CHegde, et 
al . . 19893. Figure 5.1 gives the flow path of River Ganga near 
Kanpur from 1910 to 1982. Figure 5.2 indicates the thalweg 
profile during 1947 to 1982. 

In 1910-11, a clear meander pattern was observed in 
Ganga river making a perfect northward loop between Purwa-Khargi 
and Flautapur, and another loop was observed between Rautapur and 
Bhairoghat which was reverse loop curving southwards with point of 
inflexion at Rautapur. Due to this meander pattern the main 
current up to 1912 was flowing along the right bank where Kanpur 
city is situated. Thereafter it was observed that the Ganga river 
showed signs of shifting towards the left bank. But the shift was 
very slow. Near about 194,5. the shift was more pronounced and for 
the first time the problem of water scarcity was felt at the water 
intake of city water supply work and thermal power station which 
are located on the right bank of Ganga river at Bhairoghat. 

Near about 1955 the problem was further aggravated. To 
overcome the difficulty a channel was excavated to carry a part of 
river flow from left bank to the intake of water work and thermal 
power station. This process has been continued since that time. 
But after every monsoon the channel was silted up and it has to be 
dredged. 
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Thereafter the main flow gradually shifted away from 
Kanpur and erratic movement of the meander continued. A survey 

conducted in 1982 showed a considerably transformed river 
configuration with a clearly discernible meander only upto 
Fatehpur. Downstream of Fatehpur, the river had severely eroded 
the left bank in a huge circular arc between Shankarpur and 
Gangaghat and several flow channels were formed with numerous 
shoals big and small in the bed. The transformation in the 

meander between 1910 and 1982 between Purwa-khargi to Rautapur is 
quite remarkable as would be evident from the following 
characteristics C Hegde . et al., 1989 3. 

Parameters 1910-11 1982 

Radius of curvature 7.85 Km 6.0? Km 

Meander length 12.0 Km 8.0 Km 

Meander width 1.5 Km 1.0 Km 

3.2 MORPHOLOGICAL FEATURES OF RIVER GANGA NEAR KANPUR 

The valley of the Ganga in a considerable length 
upstream and downstream of Kanpur is characterized by high firm 
bank on the right and low gently sloping bank on the left. 
Between Bithur and the railway bridge the river flows along the 
left bank. The right bank is resistant to erosion while the left 
bank, though composed of stiff hard clay, is comparatively easily 
erodible except in the reach between Rautapur and Shankarpur. 
During the last ?0 years large areas on the left bank between 
Shankarpur and railway bridge were under severe erosion at 


different times. 
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Width of the valley varies greatly, being 6.5 Km at 
iithur, 9.5 Km at Rautapur and Bhairoghat and 3.5 Km at 
arsa i yaghat . 

Downstream of Bithur, the right bank turns sharply to 
he south and the river leaves the right bank and crosses over to 
he left bank near Rautapur. This is a significant topographical 
eature that has a considerable influence on the river behaviour 
n the reach between Fatehpur and railway bridge. 

Between Rautapur and Fatehpur the river had always been 
Dre or less at the extreme left of the flood plain with the 
isult that all the spilling during flood took place on the right 
de over a large area right upto Baikunthpur on the right bank, 
lis caused building up of flood plain only on the right side, and 
e huge Baikunthpur bar gradually developed in size as well as in 
ight pushing the river more and more towards the left bank, 
ereby increasing the sinuosity of the meander. The apex of the 
ander, therefore, moved in the downstream direction. From 1953 
^ards the Kanpur side of the flood plain was gradually raised by 
rtical accretion during floods, as had been going on in the 
i kunthpur-Bha i rogha t area. This caused a fairly significant 
>pe in the flood plain towards the left bank which in turn 
ised the river to get entrenched along the left valley wall, 

oct of tho Railway Bridgo : As mentioned in the introduction 
re are sharp changes in flow direction near road and railway 
dges. These sharp changes caused excessive loss of stream 


er which 

is evident from a 

number of 

large 

shoals 

i n 

the 

i n i t y of 

the two bridges. 

In fact. 

several 

spans 

of 

the 



bridges on the left have been so badly silted up that people have 
occupied them for temporary residence. These shoals have further 
pushed the river to the left, with the result that severe erosion 
has taken place in a huge circular arc from Shankarpur ..to 
Gangaghat , 

Due to rigid control of the railway bridge in the 
downstream, there has been squeezing of the meander and thereby 
increase in the meander width between Fatehpur and the railway 
br i dge , 

The probable cause of erratic meander shift of the Ganga 
river at Kanpur according to Hedge etal. [19893 is due to <i> 
highly irregular shape of the valley in the area which caused 
development of the flood plain predominantly on one side of river; 
<ii) flood plain transformation brought about by the greater flood 
of 1924; and (iii) the location of railway and road bridges on the 
extreme right side of the flood plain. 

3.3 HYDROLOGICAL ANALYSIS OF SI REAMFLOW DATA 

For a model study of river instability and plan form 
change one needs to know the average flood and maximum flood 
hydrographs. In order to achieve this, some flow data collected 
from various sources are analysed. Ten-daily average discharge 
data of 28 years and maximum design discharge for barrage were 
available. Time series graphs have been drawn for maximum floods 
of each year and for lean period flow in each year for a period of 
28 years as shown in Figures 5.5<a> and ?.T(b). The first figure 
shows the time series of annual peak flow. From this time series 
it can be seen that maximum floods of magnitude around 9000 m^/sec 
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occurred in four years in 1961, 67, 71 and 78. But the highest 
was in 1978, around 10300 m'/sec. The second figure shows the 
time series of annual lean period flow. The lowest flow was in 
the year 1980 of magnitude around Z6.0 m^/sec. 

Ten-daily average discharge data have been used for the 
computation of the 28 year average flow hydrograph. The average 
flow hydrograph and flood flow hydrograph of 1978 have been drawn 
as shown in Figures 3.4(a> and 3.4(b>. These two hydrographs have 


been used in the model. 

The 

first 

hydrograph 

has 

been 

used 

for 

training the river, and 

the 

second 

hydrograph 

has 

been 

used 

for 


checking the stability of the river course. 

The annual peak flow data has been used to compute the 
return period, T (inverse of frequency of occurrence) of flood 
discharge. The annual flood frequency data are found to 

approximately fit a Gumbe 1 ( Extreme Value Type I ) distribution 

as shown in Figure 3.5. 

The parameters for the Gumbel distribution curve are as 
foil ows 


X = 5746,94 m"^/sec. 


Reduced variate y 


O' = 20 28.99 m /sec. 

X - 4833.81 
1582 


Cumulative probability, F<x) = exp [ - exp<-y) ] 
Exceedence probabilityP = y = 1 - F<x) 
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CHAPTER IV 


DATA ACO’lirSITTOH AND EXPERI MENI’AL DESIGN 

4^1 FIELD OBSERVATION 

Field trips were made a number of times to collect 
relevant river bed data and to note the path of flow, bed 
configurations and bed vegetations. A detailed field observation 
of the flow path, bed configurations and bed vegetation was 
carried out as first phase of investigation during the lean summer 
months . 

The first field trip was to the intake structure of 
water works at Bhairoghat situated on the right bank of Kanpur. 
The Officer Incfiarge of water intake structures informed that the 
requirement of water for water works was about 10 m^/sec per day 
and that amount of water was not available in lean seasons. This 
was due to shift of main flow current towards Unnao side. An 
unlined channel dug to convey water from the river to the intake 
structure had width about 30m at top, depth about 3m in centre, 
and length of about 6 Km. Dredging of the canal was in progress 
to keep the capacity of the canal at required level. The 
expenditure for dredging of the canal was of the order of Rs. 1.0 
lakh per month. 

Walking along the canal, one could see heaps of dredged 
sand on either side. That gave indication of the quantity of silt 
removed. Green vegetables were grown on the Kanpur side of the 
canal. Water was supplied to vegetable gardens from the canal by 
earthen pots at different places along the canal, thus damaging 



the bank and causing silting in canal. 

The second visit was to Bithur. Near Bithur the river 
divided into two streams by making a big island in between. But 
only a small quantity of water was flowing along the right bank 
and the rest of the water was flowing along the left bank. The 
farmers grew green vegetables on that island, and they had made a 
temporary bridge on the river to connect both banks. 

The third visit was to Ja 1 -Sans than , Kanpur. Hith the 
courtesy of Jal-Sansthan sufficient information such as river 
maps, topographical data and cross sectional details of the river 
Ganga were collected. 

The fourth visit was to Bithur for collection of sand 
samples. Sand samples across the width and also from the banks 
and island were collected. The locations of sand sampling at 
Bithur is shown in Figure 4.1. 

The fifth visit was again to Bhairoghat area. There is 
a big island between Bhairoghat and Unnao side bank approximately 
10 Km long and 4 Km wide. There are five villages on this island. 
Given the soil fertility the main productive activity of these 
villagers were agriculture and husbandry (milk farming). 
Villagers informed us that during high floods they all shift to 
Kanpur city temporarily and come back again after receding of the 
flood. Nine bed samples at Bhairoghat across the width were 
collected, locations of sand sampling near Bhairoghat area as 
shown in Figure 4.1. Sieve analysis was carried out to know the 
size distribution and properties of bed materials. 

The sixth trip was to the Kanpur-Lucknow road and 



Unnoo 
side bank 



SCHEMATIC OF RIVER GANGA BETWEEN BITHUR AND KANPUR 



1 ' a i 1 w ri y h r i ( J < j f * ^ i . I t w a n o b n n i' v o ci t. h a i r ) (') a i t he left side g u i d e 
bank, a big scour hole had occurred. The flow coming along the 
left bank curved steeply right skirting round the guide bank, then 
turned to tl'ie left as it approached the right bank under bridges. 
Huge amount of seoiment had deposited along the left bank side 
below the bridge. On this fertile deposition again people were 
growing green vegetables. An office of the Central Water 
Commission is situated in Champapur at Unnao side of bridge. A 
visit to the office to obtain river data was not of much use. 

4.2 BED MATERIAL DETAILS 

Sediment samples scooped from river bed was kept in 
separate polythene bags labelling them serially. At some sections 
these samples were collected below the water surface and at some 
sections on dry bar in the river bed. Siev analysis was carried 
out for each location’s samples to obtain the grain size 
d i st r i but i ons . 

Sedment sampling near Bithur : The locations of sand sampling are 
shown in Figure 4.1, Sieve analysis data are presented in Figures 
4.Z(a> to <f). It is seen from the graphs of Figure 4.2 that the 
alluvial sand is more or less uniformly graded at a spot but 
median size varies across the section. Samples 1, 5 and 6 
collected from the bed of flowing water have more or less the 
same median size. Samples 2 and 4 collected from island and 5 
from stagnated pool bed have more or less the same median size. 
The former samples which were collected under running stream have 
larger median size than those collected on the islands. 
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Figure 4-.2(a-f) Bed grain size 


DISTRIBUTION AT BiTHUR. 



For samples collected from the bed under flowing water 
the standard deviation was larger than for those samples collected 
from island. 

Sedimont sampling near Bhairoghai : Nine samples were collected, 
their locations are indicated in Figure 4.1, and sieve analysis 
data are presented in Figures 4.3<a) to <i). The graph show that 
the samples which were collected from the bed under flowing water 
have larger median sizes than those collected from the island, in 
conformity with the observation for Bithur samples. Samples 4 and 
5 have been collected from the big stable island on which 
villages are situated. The median size of these samples are 
smaller than others and with higher standard deviation for sample 
4, but lower standard deviation for sample 5- 

The grain size distribution of bed material sampled in 
the lean period indicates that material under flowing water is 
relatively coarser than that under stagnant water or on the 
exposed bars or islands. This is probably because the flood flows 
during the monsoon season which inundates the bars and stagnant 
pools carry a lot of fine, suspended sand and silt. As the flow 
subsides the fine material settles on the bed. The flow during 
the remaining period of the year is relatively free of suspended 
material, and it gradually washes out most of the fine material 
from the bed over which it flows. Thus the bed material under 
lean period flow tends to be composed of more uniform, relatively 
coarser sand, whereas the material on islands and under standing 
water tends to be relatively non-uniform and with more fine sand. 
The process may be further accentuated by fine sand, and organic 
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materials deposited by wind and overland flow. 

4*3 EXPERIMERIAL SET-UP 



A 

rectangular steel 

tank of 1.1m 

width, 

0.5 m 

depth 

and 

5 . 5m 

length 

was used for the 

exper iments . 

The 

who 1 e 

1 ength 

of 

tank 

has been divided into three parts. 

first 

and 

last parts 


having length of 25 cm and middle one with 5 * Om as shown in Figure 
4.4. The first part of tank has been used as inlet to the model. 
The last part of the modeling tank is used as sediment trap. 
Ulater was supplied from a constant level overhead tank, through a 
7.62 cm dia pipe. Water from the tank was flown into a channel 
having a U-notch fitted for discharge measurements at the end of 
the channel. A point gauge was used to record the head at the 
V-notch and another movable point gauge was used to record the 
water and bed levels on the model. In the middle 5.0m length, 
alluvial sand collected from River Ganga was filled up to height 
of 50cm as bed material of the model river. 

4.4 DETAILS OF MODEL 

It is difficult to simulate the form of the river banks, 
because of the complexity of the bank sediment properties (soil 
cohesion). Since the prototype banks are steep and stable, it was 
decided to construct banks from non-erodible material by using 
cement with sand in 1:16 proportion. For this reason, only the 
model river bed was composed of erodible alluvial sand. In the 
present study the bed was mobile for certain flow conditions. 

The length scale number Ir was fixed according to 
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the maximum space available for modeling. Length of 11 cm in 
model represented 1 Km in prototype. The vertical scale was 
obtained by measuring bed slope of model and corresponding slope 
with prototype. Same alluvial sand of the river bed has been used 
in the mode 1 . 

The tank was filled by alluvial sand upto 30 cm 
depth in the middle part and water was applied for settling the 
sand. After leveling the sand the plan form of river Ganga shown 
in Figure 3.1 for 198Z was marked on the sand. Mouldable G.I. 
sheets were lined along the bank giving steep side slope. A wet 
mixture of cement and sand in ratio 1:16 was placed behind the 
G.T. sheets. After setting of the cement sand mixture, the G.I. 
sheets were removed to provide smooth-shaped, firm banks. 
Interior details like islands and flow paths were made with 
alluvial sand, according to the 1982 survey map. 


A. 5 MODEL SCALE 


The model scales have been computed from the relations 

for mobile sediment beds. Since the same alluvial sand was used 

d 

in model river bed so grain scale numbers d^ - -g— =1. It was 
assumed that kinematic viscosity and the specific weight of water 
was same for mode and prototype. 

To maintain the same grain Reynolds number in model and 
prototype 


Re 


h d 
r r 


*r 


(I ) 
r 


1/2 


= 1 


(4.1 > 


From the grain froude criterion 
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Since Ap^ = 1, = 1, both the above criteria give the relation 


h = I 
r r 


1/2 


(4.5) 


For distorted model the Vertical scale was determined from 

h 


= S , h == I .S 
I r r r r 

r 


(4.4) 

(4.5) 

Time scale of the model is obtained from the Froude criterion 
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and discharge scale, Q = h I 
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(4.6) 


Using the above equations the following scales were obtained for 
modeling the river : 


Table 4.1 


Var i ab 1 e 

Prototype value 

Model Value 

Scale 

Length^ 

27.5 Km 

5.0 m 

9090 

S 1 ope^ 

0 . 000175 

0.0125 

0.0145 

Depth 


- 

150 

D i scharge ^ 

4800m^/sec 

0 . 0005579m^/sec 

1 . 541x10^ 

T ime 

10 days 

18.0 minutes 

799 


N.B : 1. River stretch correspond i ng to model reach. 

2. Mean river slope along valley. 

5- Maximum average flood discharge. 
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4.6 EXPERIMEMTAL PROCEDURE 

Water was fed very slowly in the tank till it matched 
the lean period flow. The discharge was increased in steps as per 
the time interval and discharge intensity indicated in the 
hydrograph. Water surface level and bed level near the guide bank 
of the railway and road bridges was measured before changing 
discharge to the next step. Thus the discharge in the model was 
varied according to the hydrograph. After attaining the average 
lean period flow, the bed profile was measured at various 
cross-sections. These bed profiles were used for drawing contours 
of bed surface, indicating the bars, islands, and other 
morphological changes. Hand feeding of sand was carried out at 
the inlet when simulating flood flows. During lean period flow 
sediment was not 'fed to the model. 

The contour of bed surface was drawn at the end of each flood 
hydrograph. The flood hydrograph used for the study, and the 
level of scoured bed and water surface recorded near guide bank at 
the end of each step of flow change, were plotted along with the 


contour maps . 



CHAPTER V 


EXPERIMENTAL STUDY OF TRAINING THE RIVER REACH 

5.1 INTRODUCTION 



aim of the 

present 

1 nvest i gat ion is 

to reduce 

scour 

around 

the guide 

bunk of 

ra i 1 way 

bridge situated 

on 

Unnao 

s i de> 

and to 

make f i otA> 

available during lean period 

near 

the 

water 


intake structures situated on the right bank near Bhairoghat. 
Both these problems can be solved by shifting the river main 
current course along the right bank from upstream of Bhairoghat to 
railway bridge. This, in fact, was the guiding principle of the 
present experiments. 

In simulating floods in the following experiments, one 
limitation faced was the lack of sediment transport data, that is 
bed load and suspended load carried by the flow. Hence as a rough 
approximation bed material was hand fed at the upstream (Bithur) 
end of the model with an attempt to balance it with the sediment 
outflow at the silt tank downstream. Average temperature of water 
was near about 50°C. 

A suntmary of the relevant parameters and the results of 
the experiments are given in Table 5.1 at this end of the chapter. 

5.2 SIMULATION OF THE EXISIING RIVER COURSE 

The model bed was prepared with islands, shoals and 
water course according to the planform given by 1982 survey map as 
indicated in Figure 5.1. The lean period flow was let first to 
test the flow course. The flood hydrograph corresponding to the ; 
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3 

average flood of 4800 m'/sec, reduced to scale is flowed in time 
steps calculated according to the time scale. Contours of the hed 
surface were drawn after measurement of bed levels during the lean 
period. The flood hydrograph and bed and water surface levels 
near guide bund are plotted in Figure 5.2(a). A photo of the flow 
path is shown in Figure 5.2(b>. The main current course was 


along the left bank from Shankarpur 

to guide 

bund . 

The 

scour at 

the guide bund 

was severe, of 

the order of 

5 cm 

in model 

(correspondingly 

6.5m in prototype). 

The 

scour 

hole 

rema i ned 


almost the same size without filling up during the lean period of 
flow. 

The maximum flood of 10240 m^/sec that occured in the 
year 1978 was simulated after remodelling the bed to the 1982 
record map. The bed levels after the lean period flow. the 
hydrograph, and water surface level and bed level records at guide 
bund are plotted as shown in Figure 5.5(a). The flow pattern over 
the whole length and at guide bund are shown in Figure 5.5(b) and 
5.5(c) respectively. The flow path adhered to the left bank from 
Shankarpur to road and railway bridge with a deep scour hole of 
6.0cm deep in model. No flow near water intake structures during 
lean period was observed. The scour hole was found to get 
partially filled with sand to near about 1.6 cm in model during 
the lean period. The flow near Bithur was divided with an island 
in between. The divided flow reunited near Rautapur. A big 
island was formed near the right bank from downstream of Bithur to 
Sarasa i yaghat near road and railway bridge. 

The design discharge of barrage 18840 m^/sec. 


was 
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allowed to be reached gradually after remodelling the bed to 1982 
records. As the flood approached the peak value, the water course 
shifted from left bank to the right bank avoiding the meander from 
Rautapur to railway road bridge. The main flow current adhered to 
the right bank making water available near at the water intake 
structure- The scour hole near guide bund also started filling up 
with sand. Now, while the main current was not hitting the bridge 
embankment, a deep pool that had been formed all along the 
embankment did not get filled up with sediment. It is possible, 
however, that this pool may get filled up by suspended sediment i'n 
prototype. Since sediment transport in model is mainly bedload, 
the process of filling up of scour hole could not be observed in 
the model. The flow course indicated in Figures 5.4(a) and 5.4(b) 
should be the path of the river to achieve the aim of its 
training. Hence permeable spurs were fixed on the river bed as 
indicated in Figure 5.5. 


5.3 TRAINING OF THE RIVER COURSE 

The flow path, to be achieved by training is indicated 
in Figure 5.4a. Permeable spurs were used as training aids to 
achieve the desired flow course. The model spurs were copper rods 
of diameter 4mm and height 150 mm. They were arranged in a row 
with lateral spacing of 4 to 5 times the diameter. It was 

observed that a single row was not effective in diverting the flow 
and reducing the magnitude of the velocity of flow which passed 
through them. In'order to reduce the velocity and to raise 
downstream bed level, a second row of spurs were contemplated. 

c : ■ -.JHY 

dix. Ho. J ■? 22 14 



Water surface 











These rods were fixed in between the spacing of the first row and 
at a distance of 7 to 8 times diameter of rod downstream of first 
row. It was observed that two rows of rods arranged in staggered 
way were effective in diverting the flow. Figures 5 - 5a and 5.5b 
indicate the two rows of permeable spurs starting from the left 
bank and penetrating into the flow at an angle. 

The hydrograph of average flood was tested on the flow 
course. After allowing the flood flow to increase and then 
decrease in time steps as computed, and on receding to the lean 
period flow the, bed levels were taken. The period of maximum flow 
was increased to test the severity of flood on the flow path and 
also on the scour at guide bund. The scour depth was reduced 
considerably. It was observed that, as flood flow increased 
the guide bund was hit by main flow current. Also, a part of • the 
flow was passing through the permeable spurs, causing a skirting 
round the guide bund. This flow pattern was considered as not 
effective. Figures 5.5 indicate the flow d i rect i ons , bed levels, 
flood hydrograph, and bed level and water surface level at guide 
bund. 

In order to achieve a better performance of the model, 
again the model was leveled according to 1982 survey map without 
removing the spurs. Again the average flood hydrograph was 
allowed to flow as indicated in Figure 5.6. While executing the 
flow of flood hydrograph and tracing the path of flow current by 
permanganate crysta 1 s , an additional row of permeable spurs near 
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Figure 5-5. one more row of permeable spurs, that is three rows of 
rods, were arranged at different lengths as indicated in Figure 
5.6. The scour was reduced considerably at guide bund and the 
performance of this arrangement was found to be satisfactory. 
However, it was noticed that the length of the first array of 
permeable spur was rather excessive. 

Aiming for further reduction in the installation of 
spurs, the model indicated in Figure 5.7 was conceived in 
order to divert the water course right from Bithur. With this 
consideration the spurs were arranged in the way shown in 
Figures 5.7(a) and 5.7(b). The simulated flood was that of 1978 
(10240 m^/sec). The bed levels, flow paths, and scour details 
near guide bund were plotted. The flow path was along the right 
bank from upstream of Bhairoghat to road and railway bridges. 
There was a remarkable reduction in scour depth near the guide 
bund. Further, along the whole of left bank from upstream of 
Rautapur to road and railway bridges, a bar developed due to 
sediment deposition. This part of the island in prototype may be 
reclaimed for use. Along the right bank only one bar formation 
from downstream of Bithur to far upstream of Bhairoghat was 
observed. The river flow path from Bithur to the road and railway 
bridges was almost straight during flood, and meandering during 
lean period. Meandering width was also reduced considerably. 
Another interesting result of the above run was that the flow was 
brought near jthe water intake structure during the lean period. 



Water surface 



Figure 5.6 River training using permeable spurs Q 
4-830 M^/SEC “ 0 00036 m^/sec). 




array was removed. The bed was modelled according to the 1982 
survey map without removing the other spur arrays, and subjected 
to average annual flood hydrograph. The bed level. flow path, 
islands, shoals and scour details near guide bank and the flood 
hydrograph are shown in Figure 5.8<a). Figure ^.8<b) indicates 
the flow pattern during transition of shifting from left bank to 
right bank near Bhairoghat. The flow path during lean period was 
at the water intake structures, and all along the right bank to 
bridge passage. The scour depth was reduced considerably. The 
earthen embankment built to protect the road and railway bridge 
was thus completely safe from the flood attack. Further, a big 
bar on the left bank was formed, indicating the possibility of 
reclaiming this land. This configuration of permeable spurs 

arrangement was now simulated with the 1978 flood hydrograph. 

5.4 CHECKING FOR THE 1078 FLOOD 

The flood hydrograph corresponding to the 1978 flood was 
allowed to flow through the model river with spur arrays as 
indicated in Figure 5.8. After the run the bed levels, flow paths 
during lean period, islands, hydrograph, and bed and water surface 
levels at bridge guide bank were plotted as indicated in Figure 
5.9. Figure 5.9 indicates that the performance of permeable spur 
was better in this case than for the average annual flood. The 
scour depth was reduced considerably, the island width on the left 
bank increased and the lean period flow was along the right bank, 
from far upstream of Bhairoghat till the road and railway bridges. 



Water surface 



Figure 5.8Ca) River training using permeable spurs, Q 
4-830 M^/sEC “ 0-00036 m^/sec). 
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From the above test it was concluded that the permeable 
spur arrangement indicated in Figure 5.8 is adequate for training 
the river to achieve two objectives, namely to divert the lean 
period flow near the water intake structure and to reduce the 
current attack and scour at the guide bank provided to protect the 
road and railway lines at Unnao side. Table 5.1 summarizes the 
comparative features of different experiments conducted in this 


study . 



Tablo S- 1 Sammairy Information of Models Tested 


Mode 1 

Flood 

Train i ng 


Model Results 

run 

hydrograph 

works 

scour at 
embankment 

planform and other 
changes 

1 . 

Average 
annua 1 
flood , 

4820 m'/sec 

None 

5.0 cm at 
peak flood 
and reduced 
by 1.0 cm 
in lean 
period. 

Some islands washed out 
between Rautapur and 
Fatehpur and formed 
small bars. Flow path 
was along the left bank. 

2 

1978 flood 

3 

= 10240 m /sec 

None 

6 . 0 cm at 
high flood 
and reduced 
by 1.6 cm 
in lean 
period. 

Flow path divided by 
island at Bithur. Many 
shoals and bars between 
Fatehpur and Shankarpur 
and a big island formed 
near right bank downstrea 


of B i thur . 

Flow path comes near wate 
works in high flood but 
along left bank in lean 
period. 


Des i gn 

None 

5 . 0 cm and 

Big island scoured from 

d i scharge 


reduced by 

right bank and bars 

of barrage 


2.5 cm 

formed along left hank. 

= 18840 m^/sec 


in lean 
per i od . 

Flow shifted to right 
bank in lean period 
also. 

Average 

Doub 1 e 

2 . 0 cm and 

Flow was divided by 

annua 1 

row of 

compl et e 1 y 

island at Bithur and 

f 1 ood , 

spurs 

reduced 

big bars formed behind 

= 4820 m^/sec 

4 to 5 

in lean 

spurs . 


dia c/c 

per i od , 

Flow along right bank 


near left 


but direct hit to guide 


bank close 

bank and embankment at 


to Rautapur 

1 ef t s i de . 

Average 

Three 

1.5 cm and 

Big sand bars formed 

annua 1 

rows of 

compl ete 1 y 

along the left bank 

flood , 

spurs a 

reduced in 

from Rautapur to. 

= 4820 m'^/sec 

4 to 5 

1 ean 

Shankarpur; small sand 


dia c/c 

period. 

bars along right bank. 


wi th 


Flow along right bank 


addi t i onal 


and some flow divert 


spurs at an 

towords left in lean 


in front of 

period. 


guide bank 

• 

: j 
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6 . 

1978 

Double 

2.0 cm and 

Flow divided at Bithur 


flood , 

row of 

completely 

and island formed down- 


= 10240m^/sec 

spurs 

reduced in 

stream of Bithur. Big 



a 4 to 

lean period 

sand bars along left 



5 d i a 


bank from upstream of 



c/c near 


Rautapur to upstream 



left bank 

■ 

of railway bridge. 

7. 

Average 

Doub 1 e 

^.0 cm and 

Flow path was along the 


annua 1 

row of 

reduced by 

Right bank during lean 


flood , 

spurs a 

2.0 cm in 

period. 


= 4820 m^/sec 

4 to 5 

lean period 

Flow divided by island 



dia c/c 


at Bithur; big sand 



near left 


bar formed along left 



bank 


bank, also along right 



and along 


bank . 



the guide 


Flow completely shifted 



bank . 


toward right bank in 





lean period also. 

8 . 

1978 

-do- 

1.0 cm and 

Sand bars formed along 


flood 


comp 1 e te 1 y 

whole left bank from 


=10240 m^/sec 


filled up 

upstream of Rautapur to 




in lean 

upstream of railway 




period 

bridge. Bars also 


existing at right 
bank . 

Flow along the right 
bank making water 
available at water 
works, round the 
year . 



CHAPTER VI 

CONCLUSIONS, LIMITATIONS AND RECOMMENDATIONS 


.1 CONCLUSIONS 

Hydraulic modelling of River Ganga near Kanpur was 
;arried out to achieve two objectives namely <a> to bring the flow 
luring surrmer months near Bhairoghat where city water supply 
intake structure is situated, and (b) to reduce the severity of 
scour occur ing at the bend near Champapur, which has an earthen 
embankment built to protect the railway and road bridge connection 
Kanpur to Lucknow. After full reconnaissance of the river reach 
near Kanpur and collection of flood data and bed level data from 
different sources, a distorted hydraulic model with horizontal 
scale 1 in 9090 was prepared. Tfie average annual flood and the 
peak annual flood of 28 years were used for simulating the river 
discharge in the model. The model was tested for its ability to 
reproduce the existing plan form pattern of the river using the 
average annual flow hydrograph. The highest flood intensity used 
for the design of the proposed barrage by Kanpur Development 
Authority, indicated the flow path for which the river has to be 
trained. 

Permeable spurs made of copper rods of diameter 4 mm 
spaced laterally 4 to > diameter apart in two rows in staggered, 
arrangement, with spacing between the rows of 8 diameter was 
adopted to train the flow path in the model. After various 
trials, the spur arrangement indicated in Figure 9.8 was evolved 
to achieve the two aims cited above. This spur arrangement was 
tested for the highest flood in 28 years, wherein its performance 
improved in comparison with the 28 years average flood hydrograph. 



T both cases the lean period flow skirted the water intake works 


t the river’s right bank, and considerably checked water velocity 
id scour near the embankment. The river thalweg had a stable 
rofile closer to the right bank. 


.2 LIMITATIONS 

The main limitations of the present model are the 
dl lowi ng ; 

(a) The horizontal scale of the model is too large and its 
vertical scale is also large. The ratio of horizontal 
to vertical scale should not exceed about 6:1 for 
reliable results. 

(b) Materials like coaldust, which are lighter in density 
in comparison to alluvial sand, should be used in the 
model bed to simulate bed mobility with less scale 
distortion. 

<c> The sediment transport in model is by bed load only 
whereas in prototype, the sediment is moved mostly in 
suspension. This can be achieved by changing to less 
dense bed material. 

<d> The sediment load data of the river were not available. 
Efforts should be made to get these data from Central 
Water Commission Office situated near Champapur or 
from other sources. 


6.3 RECOMMENDATIONS FOR FURTHER WORK* 

A hydraulic model to a reasonable horizontal scale of 1 
in 600 and verti,iai scale of lin 100 with coal dust as the bed 



iterial should be built for detailed study. An 
jpraisal of the use of permeable spurs as river training 
)mparison to building the proposed barrage at Bhairoghat 
j made for appropriate selection of the one that best su 


economi c 
aid in 
should 
its the 


■oi n ob j ect i ves . 
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APPENDIX A 

A.l DIMENSIONAL ANALYSIS OF SEDIMENT MOTION : 


For studying the very complicated processes at the initiation 
of sediment motion and during sediment transport, where the 
transporting fluid and the transported sediments interact with 
each other, it is advantageous to consider the various factors of 
influence beforehand. 

The specific sediment transport rate . expressed as mass 
per unit time and unit width, depends upon the material properties 
of the fluid and the sediment, upon the grain diameter d of the 
transported material as well as upon the water depth h and the i 
slope S of the transporting flow. 

f(g . p , , g , p . d . h . S ) = 0 (1) 

By dimensional analysis <n theorem of Vaschy), five 
dimensionless parameters can be derived from this equation. The 
sediment transport is accordingly characterized by 


= f 


Fr. 


Re 


w 


s w 


h 

d 


< 2 ) 


In this 
un it t i me i s 
number 


equation, the sediment transport per unit width and 
expressed by the dimensionless sediment transport 



< 5 > 


As a reference velocity the shear velocity as defined by 

Prandtl is used, which is defined by the bottom shear stress and 



the density of the fluid. 




» - [r^/p = jifT 


(4) 


The parameters 


Fr 


w 


y -y 
^ s w 


he. w 

3 ' 

d 

s w 


V 

* 

gd 


Fr 


* (y -y 
s w 


<5> 


is a modified grain Froude number, containing the ratio of the 
specific weights of the fluid and the submerged sediment. It 
corresponds to the shear stress parameter as defined by Shields. 


The Reynolds number of the grain 


Re, 


■4 ghs d _ 


w 


w 


(6) 


has been introduced by Shields in 1956 in this form. 


By connecting the modified grain Froude number with the 
Reynolds number of the grain, a further dimensionless quantity can 
be defined. 

2 

Dividing Fr^ by Re^ 


Fr. 


Re' 


w 


V 


<y -y ) 
s ' w 


w 


gd' 


(7) 


The expression can be considered as a dimensionless buoyancy 
parameter, it contains exclusively the properties of the fluid and 
of the sediment and characterize the material properties of both 


media. 


t 


A. 2 SEDIMENT TRANSPORT EQUATIONS : 


Shields formulated from his experimental values the following 
relationship for sediment transport 


q (j/ > 

s s w 
q r s 


10<T -T ) 
s c 

-p >• d 


TP 


s 


w 


< 8 > 


in which and q represents the sediment and water volume per 

unit time and unit width. 

Above equation can easily be transferred into the 

dimensionless quantities defined above. 




^ 10 Fr^(Fr^-Fr^^> 
* 


< 9 > 


the term <V /V„) represents the frictional resistance of the bed. 
This resistance depends upon Re^ and the ratio of the water depth 
to the roughness elevation or grain size <h/K^) or (h/d) . 


H. A. 

Einstein 

i nt roduced 

the 

concept 

of 

probab i 1 i ty 

for 

anal yz i ng 

sed iment 

motion. On 

the 

bas i s 

of 

observat i ons 

and 


probabilistic considerations, he found that the sediment transport 
depends upon the "transport intensity". 



( 10 ) 


and the "flow intensity" 


V' 


P ~P 

w 


w 


C Y. ) S 


’y s 


<11 ) 



The relationship between the parameters (p and ^ is given in 
f i gure < A . 1 ) . 


Einstein bed load function can easily be expressed in terms 
of the dimensionless parameters defined above. 


For the transport intensity (p , 



( 12 ) 


Reformulation of the flow intensity yields 


V' 


iy -y y d, 
s w 



( 1 ?) 


with these expressions, the Einstein bed load function can also be 
presented as a function of g^ and F r^ . 

A. 3 DERIVATION OF MODEL SCALES 


It has been seen that the sediment transport is dependent 
upon the dimensional parameters Re^ and Fr^. 


If sediment motion is to be modelled correctly in a hydraulic 
model, then the model sediment and the model scales have to be 
chosen such that these dimensionless parameters have the same 
value in model and in nature. Furthermore, it is assumed that 
Froude similarity exists with regard to the hydraulic processes. 
If hydraulic similarity according to Froude’s law is given, then 
one can derive the scale relationships by equating the Reynolds 
number of the grain and the modified grain Froude number in nature 


and mode 1 . 
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(14) 
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Re 
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Re 
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Fr 
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Fr 
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For the derivation of model scales, one has to introduce now in 
addition to the length scales for horizontal lengths and h^ for 
vertical lengths two additional spale numbers. The grain size 
scale number 

d = d /d (1^> 
r n m 


and the scale number of the sediment density (under submerged 
cond i t i ons > 


Ap 




n 




m 


ip -p > 

w n 

ip -p ) 

*^3 w m 


iy -y ) 
s ^ w n 

) 

s w m 


(16) 


If the kinematic viscosity and the specific weight of water 
is assumed to be equal in model and nature, then one obtains by 
equating the Reynolds numbers of the grain for model and prototype 


Re 




g 1/2 h 1/2 S 1/2 d 
r r r r 


x> 
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h 1/2 d 
r r 

h '^2 
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h d 
r r 

TTTT 


(17) 


Similarly, the modified Froude number in model and prototype can 
be equated if the specific weight of water is the same in model 
and nature, then one obtains with equations 


Fr 


*r 


p g 
wr r 

Ap g 
r r 
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Kp d n 

r r 
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r r 



(18) 


The scale relationships equations (17) and (18) yield in 
combi nat i on 


h 
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5 


Ap 

r 


Z 

r 


5/2 


(19) 



This equation resembles the scale relationship 
by Lacey and Inglis for models with mobile beds 


found empirically 



In order to determine the correspondence between the model 
scales and the roughness related, to the grain diameter, the 
Mann i ng-Str i ck 1 er equation can be used. For the roughness 
coefficient according to Mann i ng-St r i ck 1 er , 


1/2 


<K 


^ r 


<h ) 
r 


2/5 


( 20 ) 


For a mobile bed, the roughness coefficient can be 

determined from the grain diameter d^^ (90^ sieve passage) of the 
cover layer according to the empirical relation. 


= 26/dgQ^^^ and with that 


Now simple relationship can be derived between the horizontal and 
the vertical length scale 



( 22 ) 


A relationship between the scale number for the density of 
the model sediment and the scale of the grain diameter can be 
der i ved 

d = Ap (25) 
r r 


The time scale for the sediment transport can be determined 



Dy equating the dimensionless sediment transport number in 
nature and model. If the transported sediment rate . expressed 
as dry weight per unit width and unit time, is divided by then 
one obtains the corresponding volume per unit width and unit time 
without void space. For sedimentation and erosion processes. one 
has therefore consider in addition to the natural structure of the 
bed, which depend's on the grain sieve curve and other factors by 
introducing the void fraction p. relationship results 


* r r 


4 h 

f r r 


. , 1/2 . 1/2 c 1/2 

t . d . g h .5 

s , r r r r r 


= 1 


and therefore t 


s , r 


5/2 


“ Pr 


(24) 


For identical sieve curves in nature and model one can 


approximately assume <p^ = 1). 


Then 


s , r 


. ^/2 A 

4r 

TTY” 


<25) 


If the flow is assumed to follow Froude’s law model law, then the 
time scale for hydraulic processes 


t 


h » r 



(26) 


The time scales for sediment transport processes and for 
hydraulic processes in a model are not the same. 


Av«r«g# T#n daily Di«charg# in cum^cm 


Riv#r GANCA 


Sit# KANPUR 




I - 

1 to 10 

. II - 10 to 20 

, in- 

20 to Last 


Hont h 

1 960 



1961 



1962 



I 

II 

III 

I 

II 

III 

I 

II 

III 

Jan 

153.6 

148.39 

151.4 

269.42 

316.05 

249.24 

374.56 

314.18 

317.95 

Feb 

159.7S 

142.56 

114.27 

384. 04 

707.52 

458.94 

565.55 

420.94 

328.01 ^ 

March 

93 . 8 

129.7 

165.1 

266,43 

202.65 

167.26 

303.77 

537.46 

397.52 

Apr i 1 

125,88 

91.38 

80.85 

148,24 

122.17 

105.48 

247.47 

176.26 

132.52 

May 

60.5 

52. 1 

42. 12 

92.95 

81 .69 

72.73 

151.19 

135,35 

96.33 

June 

43 . 93 

47.43 

180.38 

64.73 

394.98 

633.53 

93.7 

141 . 62 

397.44 

July 

422 ,81 

3956 . 67 

5281 . 15 

759 . 64 

371 0 . 86 

3136.5 

678.88 

908.72 

2299. 1 

Aug 

3671 . 1 6 

4812.76 

4641 .81 

6732,39 

7087.24 

9144.19 

4576.85 

3916.92 

4040 . 13 

Sept 

4327. 14 

5238.75 

1959.79 

5974. 16 

3147.31 

8069.46 

4231 .56 

2581 . 02 

4902.62 

Oct 

5355.63 

4120.64 

1128.87 

1874.38 

4368.41 

1805.73 

2187.49 

978.49 

666.84 

Nov 

863.8 

695.34 

448. 06 

1271 .24 

922,85 

767.31 

560.83 

466.62 

313.5 

Dec 

321.91 

287.77 

241 . 99 

551 ,86 

392.68 

429. 17 

252.46 

240.86 

220.75 


Month 

1 963 



1964 



1965 



I 

II 

III 

I 

II 

III 

I 

II 

III 

Jan 

200.69 

184.61 

230 . 04 

203.07 

173.74 

175-34 

211.41 

174,13 

178.35 

F eb 

197.93 

154.02 

130.55 

160.50 

155,18 

140.89 

157.17 

141.6 

162.61 

March 

131 .72 

129.95 

132.81 

136.82 

120.63 

135.83 

153.42 

127.87 

149.39 

Apr i 1 

135,67 

1 07.41 

90 . 7 

140,48 

103.02 

79.82 

194.03 

174.17 

97.21 

May 

70,3 

72.24 

55. 8t 

71 .92 

74.98 

58.22 

80,68 

64.19 

63.6 

June 

55.21 

276 , 1 7 

501 .54 

53.2 

49.2 

63.31 

55.9 

56 . 04 

1 03.08 

July 

627.57 

1288.75 

1823. 14 

469 . 06 

3182.98 

3404.93 

411.03 

666.61 

1849. 03 

Aug 

3127.01 

4473.51 

7764 . 92 

3684.48 

3473.97 

3789.24 

2930,46 

1 983.66 

1 600 . 95 

Sept 

6987 . 03 

4886.2 

6703 . 02 

4591 . 99 

4285. 81 

3014 . 95 

2277,91 

2396.4 

923.97 

Oct 

1901.07 

942.89 

696.35 

3295.05 

1262. 64 

936.07 

430.2 

319.63 

£70. 14 

Nov 

673.99 

714.76 

439.29 

760.77 

581 . 94 

443.67 

181 .31 

167. 07 

151 .29 

Dec 

291 .43 

260. 1 3 

240.38 

307.49 

326 . 46 

290.44 

129.06 

1 08.63 

109.17 


Month 

1 966 



1967 



1968 



I 

II 

III 

I 

II 

III 

I 

II 

III 

Jan 

95.48 

74.4 

75 . 02 

114.16 

111.23 

94.23 

373.39 

373.44 

321.74 

Feb 

81 . 94 

81 . 89 

92.71 

85.21 

84.18 

73.38 

266.66 

226.85 

200.79 

March 

83 . 84 

67 . 27 

55 , 08 

67 . 64 

68.33 

68.64 

220 . 43 

174,59 

168.16 

Apr i 1 

46 . 0 

43.48 

39 .69 

79 . 43 

62.26 

54.59 

181,94 

122.75 

139.48 

May 

36 . 69 

43.3 

56.26 

47.94 

42.84 

40.39 

151.19 

125.97 

78,48 

June 

34.18 

29 . 76 

301 .58 

28.33 

173.01 

314.23 

56.75 

843.31 

805.51 

July 

2392.84 

821 . 38 

1 955 . 05 

078 . 25 

1874 .2 

2931 . 73 

1848. 14 

2271 . 09 

4247.46 

Aug 

4875.73 

7100.54 

4559.86 

3344 . 83 

6099 .57 

6775. 13 

2622.55 

3908.74 

4077.30 

Sept 

2301 .4 

t 951 . 07 

1254.27 

9232.82 

5338.9 

2524,33 

1681 .91 

1 120 .29 

1972.63 

Oct 

695.64 

483.84 

317.43 

1257.25 

726 .18 

444.2 

887.74 

498.56 

351 ,23 

Nov 

204.38 

180.18 

147.31 

352.66 

421 .2 

348.65 

265.14 

233. 16 

203.89 

Dec 

141.32 

136 . 92 

123.06 

267.82 

256.57 

250.22 

1 69,87 

136 , 06 

124,38 



Month 


1969 1970 1971 



I 

I I 

I I I 

I 

1 1 

III 

I 

1 1 

III 

Jan 

180.54 

122.81 

123.65 

171.81 

153.04 

195.13 

122.32 

116.06 

116.33 

Feb 

116.53 

114.51 

too. 6 

236.54 

197,71 

221 . 73 

118.23 

101.92 

93.68 

March 

91 .51 

82 . 44 

8 0.41 

279.59 

297.52 

230.8 

1 00.48 

111.29 

89.61 

April 

69.33 

59.14 

56.02 

175. 02 

134.3 

115.24 

75.74 

68.04 

80 . 04 

May 

52.76 

57.83 

45.27 

95.35 

64.23 

65.56 

68.17 

75.66 

74.73 

June 

54.38 

271 . 97 

293.59 

86.24 

341 .85 

621 . 96 

148.94 

1187. 07 

1770.31 

July 

468.39 

1292.79 

2573. 18 

960 . 72 

1583.87 

1797.97 

3992.43 

3772.81 

3065.81 

Aug 

3263 . 95 

3883. 61 

5176. 1 

2863.26 

4830 42 

3168.51 

5481 .77 

8863.49 

5170.5 

Sept 

2454 . 06 

3367.64 

5388. 17 

2912.77 

3115.72 

1917.63 

71 03.0 

5630.05 

1718.86 

Oct 

3637.56 

1115.8 

755. 64 

1315.37 

987.21 

666 . 63 

1147.65 

826, 16 

1966.06 

Nov 

828.86 

562.63 

386,73 

451 . 74 

281.5 

232.94 

1015.81 

748.89 

548.01 

Dec 

316.84 

260.45 

206.28 

199.03 

167.78 

144.17 

393.87 

536.04 

£85.52 


Mont h 

1 972 



1 973 



1974 



I 

II 

III 

I 

II 

III 

I 

II 

III 

J an 

241 . 68 

232. 19 

208.06 

141.73 

137.07 

139.15 

205,23 

100.72 

157.73 

Feb 

asB . 1 7 

307.53 

264.49 

149.93 

122.78 

106.89 

135.59 

127.47 

122.79 

March 

t 95.29 

171.17 

1 31 .76 

10 6.18, 

94-71 

100.0 

1 07.51 

99.64 

95.3 

Apri 1 

121 ,99 

1 07 38 

1 02.59 

80.36 

79 . 24 

64.62 

100,7 

101.34 

80.61 

May 

105.6 

76 . 97 

67 . 04 

69 . 9 

208.6 

128.92 

69.91 

115.45 

95.47 

June 

103.65 

111.72 

810.01 

128.06 

449.95 

2798-67 

67.61 

59,99 

157.39 

July 

753.44 

2906 . 12 

1820 , 14 

2079.43 

2248. 14 

4245.06 

202.52 

490.45 

2490.94 

Aug 

1302.93 

1918.32 

2647.32 

6688.8 

6742.65 

5000 . 64 

3781 . 01 

4287.67 

1795.38 

sept 

2540.53 

5530 . 49 

3399 . 76 

3426 . 39 

3235.87 

2197.3 

1684-94 

855.22 

573.88 

Oct 

1 153.53 

687.77 

497.28 

2068,54 

1933.64 

1 004.28 

444- 08 

283.82 

839 . 12 

Nov 

477.36 

321 , 27 

251 . 66 

652.42 

468.63 

367.08 

162.23 

132.5 

116.29 

Dec 

230.84 

21 6 . 51 

1 99 . 93 

310,79 

272.61 

231 .73 

100.96 

90.07 

97.2 


Mont h 

1 975 



1976 



1 977 



I 

n 

III 

I 

II 

III 

I 

II 

III 

Jan 

95.13 

87.24 

79 . 05 

168.07 

154.67 

141 .64 

126.4 

108.06 

121 .56 

Feb 

70 ,31 

68.67 

67.12 

136.45 

122.23 

137.93 

156.45 

146.74 

1 83 . 43 

March 

63.39 

55.86 

51 .38 

131 .8 

102.58 

97.81 

116.27 

96.47 

85. 13 

Apr i 1 

46.74 

43.13 

40.14 

91.91 

84.19 

89.07 

77.28 

80.35 

69. 52 

May 

33 . 46 

40 . 64 

91 . 9 

106.2 

77.52 

100.54 

68.19 

64.3 

63,3 

June 

92.11 

129.55 

1284.25 

143.14 

464.6 

476 . 1 6 

62.83 

54.23 

67. 15 

July 

1476 .4 

1423,38 

3106.9 

420 . 66 

1643.94 

2694.75 

515.27 

3359.22 

4796.81 

Aug 

4979 .42 

3529.32 

3577.31 

3607. 1 t 

3899.34 

5974.19 

4872.58 

4126.9 

8972.85 

sept 

3483 . 08 

6705.69 

3561 .91 

2812.14 

2233.79 

1 191 .87 

3555.53 

4024.83 

31 17.55 

Oct 

1790 , 88 

1059.43 

760 .73 

783.02 

500.37 

289.32 

1704.74 

1 024.71 

664.63 

Nov 

564.87 

350.81 

283.27 

353.24 

297.76 

195.61 

445.38 

351 .5 

295.8 

Dec 

239 . 93 

215.15 

1 94 , 78 

164.92 

154.92 

146.47 

256.28 

216.76 

195.83 


Mont hi ' ’ 

’ 1-978 



1979 



1 98 0 



I 

II 

III 

I 

n 

III 

I 

II 

III 

Jan 

205 .43 

183.88 

195.0 

835.91 

816.15 

353.48 

77.43 81,7 

88.78 

fmb 

IBS. as 

174.38 

841 .44 

317.77 

348.08 

397.8 

71 .95 

: 65.86 

68.58 

Harch 197.86 

884.76 

4 0 6.90 

407.58 

324 . I 8 

869.75 

56.87 

59.14 

50. 1 8 

Apr i 1 

ai8 .5 

147.79 

148.96 

178.85 

1 68 . 8 

874.55 

58 . 93 

51 .48 

54 . 9 

Hay 

1 30 . 09 

144.83 

844.6 

880.3 

848 .18 

174.81 

41 .87 

87.59 

86.85 

Jun# 

309.44 

488,74 

1 588 . a 

133.51 

185.18 

316.83 

30.8 

133.48 

353.65 

July 

3Sa7.77 

3968. 04 

3905.91 

t 015, 78 

1363.76 

8416.3 

1 087.3 

3058.99 

5359.31 

Aug 

5576. a ' 

10838.31 

6470.86 

8159.04 

8381 .47 

8799.4 

6333.86 

5434. 63 

8485,08 

S#pt 

7309 . 48 

5469.67 

2734.42 

1059.0 

650.8 

373,01 

8858.05 

1860.57 

1385.85 

Oct 

1506.44 

805. 06 

774. 15 

814.47 

140.97 

116.04 

753, 16 

494.88 

398.81 

Nov 

689.41 

■ 489 . 54 

380.54 

97.35 

90.43 

98.38 

390.9 

838. 14 

180, t 5 

Dec 

319.14 

385.89 

313.9 

97.9 

89.89 

85.83 

166.03 150.89 

152.69 

Hont h 

1 981 



1988 



1 983 



I 

II 

III 

I 

II 

III 

I 

n 

III 

Jan 

155.84 

135.83 

130.19 

106.59 

189,30 

189.65 

160. 58 

157. 16 

134,85 

Feb 

148.36 

181.86 

102.88 

160.56 

137.68 

188.85 

807.54 

885.50 

143.68 

Harch 

89,46 

96.87 

106,39 

114.46 

1 96.81 

166.87 

113.46 

96,86 

96.27 

April 

03.55 

108.33 

118.88 

163.80 

119.70 

106.33 

88.88 

86. 18 

277.64 

Hay 

71 .83 

65.90 

6 0,14 

805.96 

835.51 

877.74 

851.11 

179.67 

839,71 

June 

88.74 

108.90 

99.83 

153.95 

1 89,16 

499.01 

301 . 05 

308. 05 

867,26 

July 

1491 . 08 

1 704 . 74 

3587.59 

514.42 

640.38 

1 978.88 

1084.35 

936.43 

1130.77 

Aug 

5889 . 78 

3374. 09 

8489.85 

4141.89 

3930 . 97 

5673. 15 

8667.38 

3777.59 

4531.88 

Sept 

1480.47 

771 . 00 

554.85 

6484.88 

3889.88 

1804.38 

4508.89 

5665.68 

5569.90 

Oct 

663.94 

398.68 

880.84 

560.69 

319.64 

846,35 

3401 .89 

1946.88 

984.64 

Nov 

196.38 

880. 13 

197.51 

838.70 

819 .83 

198.41 

647.08 

438. 16 

336.99 

Dec 

155.34 

135.70 

117.80 

178.44 

138.10 

153.70 

887.66 

870.38 

260.78 

Hont h 

1 984 



1985 



1986 



I 

II 

III 

I 

II 

III 

I 

II 


Jan 

874.41 

885,84 

840.67 

187.41 

123.67 

123.86 

369.44 

885 .17 

. 823,81 ; 

Feb 

191,87 

186 . 90 

888.80 

115.01 

100.88 

94.08 

839.44 

331 .70 

361.99 

March 

818.64 

143.95 

185.19 

83.31 

74.80 

71.10 

876 . 99 

1 97.38 

106, 40 

Apr i 1 

138,43 

184.98 

185.61 

74.13 

64.44 

63.60 

171.19 

124.14 

104.60 

May 

101.07 

133.50 

116.45 

60.74 

59.14 

58.85 

108.77 

91 .67 

141.80 

June 

158.06 

599 . 08 

1014.08 

60 .01 

65. 16 

67.58 

113.19 

93.34 

467. 91 

July 

1548.78 

1780 .89 

8818. 65 

157.76 

586 . 69 

8317.18 

809.68 

88 18. 07 

4946, 00 

Aug 

3080.30 

1559.32 

8898.94 

8480 . 18 

8939.80 

5096,97 

4458.70 

3606.90 

4957.01 

Sept 

4853. 18 

8755.87 

1895 . 08 

3889 . 94 

3589 . 56 

4016.93 

1 994.67 

1 151 .80 

784.89 

Oct 

649.88 

541 .80 

430 . 74 

1786 .87 

5078,88 

3136.65 

1153.41 

718.48 

569.58 f 

Mo V 

855.98 

197.31 

167.40 

1864.55 

768. 14 

541 .84 

438.50 

308.58 

889.07 i 

Dec 

1 65.87 

143.86 

138.93 

408.78 

887.86 

375.17 

179.50 

813.86 

818.78 : 




t 987 





i 

f 


I 

II 

III 


I 

II 

in 


E 

[ 

Jan 

157.98 

150.79 

166,68 

July 

1 08.84 

400.49 585. 

59 


Feb 

158.31 

137.98 

147.36 

Aug 

901 .08 

744.79 708. 

71 


March 

148,44 

185.45 

111.88 

Sept 

1884.40 

1386.93 987. 

79 


Apr i 1 

1 08.63 

87. 1 

71 .34 

Oct 

643.16 

183.08 140. 

66 

■ ' ! 

Hay 

88.05 

97.14 

78.18 

Mov 

1 03.68 

118.08 91. 

64 

j; 

J u n e 

64.94 

146.55 

1 85 . 55 

Dec 

76.55 

71.66 68. 

81 




